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Abstract 

We show that, for the purpose of pricing Swaptions, the Swap rate and the corresponding Forward rates can be 
considered lognormal under a single martingale measure. Swaptions can then be priced as options on a basket of 
lognormal assets and an approximation formula is derived for such options. This formula is centered around a Black- 
Scholes price with an appropriate volatility, plus a correction term that can be interpreted as the expected tracking 
error. The calibration problem can then be solved very efficiently using semidefinite programming. 
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1 Introduction 

In the original Black & Scholes (1973) model, there is a one-to-one coiTespondence between the price of an option 
and the volatility of the underlying asset. In fact, options are most often directly quoted in terms of their Black & 
Scholes (1973) implied volatility. In the case of options on multiple assets such as basket options, that one-to-one 
correspondence between market prices and covariance is lost. The market quotes basket options in terms of their 
Black & Scholes (1973) volatility but has no direct way of describing the link between this volatility and that of the 
individual assets in the basket. Today, this is not yet critically important in equity markets where most of the trading 
in basket options is concentrated among a few index options, we will see however that it is crucial in interest rate 
derivative markets where most of the volatility information is contained in a rather diverse set of basket options. 

Indeed, a large part of the liquidity in interest rate option markets is concentrated in European Caps and Swaptions 
and, as always, market operators are faced with a modelling dilemma: on one hand, the arbitrage-free price derived 
from a dynamic hedging strategy a la Black & Scholes (1973) and Merton (1973) has become a central reference in 
the pricing and risk-management of financial derivatives, on the other hand however, every market operator knows 
that the data they calibrate on is not arbitrage free because of market imperfections. Beyond these discrepancies in the 
data, daily model recalibration and the non-convexity of most current calibration methods only add further instability 
to the derivative pricing, hedging and risk-management process by exposing these computations to purely numerical 
noise. One of the crucial filters standing between those two sets of prices (market data and computed derivative prices) 
is the model calibration algorithm. 

Recent developments in interest rates modelling have led to a form of technological asymmetry on this topic. The- 
oretically, models such as the Libor market model of interest rates (see Brace, Gatarek & Musiela (1997), Miltersen, 
Sandmann & Sondermann (1995)) or the affine Gaussian models (see El Karoui & Lacoste (1992) or Duffie & Kan 
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(1996) ) allow a very rich modelling and pricing of the basic interest rate options (Caps and Swaptions) at-the-money. 
However, due to the inefficiency and instability of the calibration procedure, only a small part of the market covariance 
information that could be accounted for in the model is actually exploited. To be precise, the most common calibration 
techniques (see for example Longstaff, Santa-Clara & Schwartz (2000)) perform a completely implicit fit on the Caplet 
variances while only a partial fit is made on the correlation information available in Swaptions. Because of these Um- 
itations, a statistical estimate must be substituted to the market information on the Forward Libors correlation matrix 
as the numerical complexity and instability of the calibration process makes it impossible to calibrate a full market 
covariance matrix. As a direct consequence, these cahbration algorithms fail in one of their primary mission: they are 
very poor market risk visualization tools. The Forward rates covariance matrix plays an increasingly important role in 
exotic interest rate derivatives modelling and there is a need for a calibration algorithm that allows the retrieval of a 
maximum amount of covariance information from the market. 

In the Libor market model, we write Swaps as baskets of Forwards. As already observed by Rebonato (1998) 
among others, the weights in this decomposition are empirically very stable. In section two, we show that this key 
empirical fact is indeed accurately reproduced by the model. We then show that the drift term coming from the 
change of measure between the forward and the swap martingale measures can be neglected in the computation of the 
Swaption price, thus allowing these options to be priced using the lognormal approximations first detailed in Huynh 
(1994) and Musiela & Rutkowski (1997). In particular, this will allow us to reduce the problem of pricing Swaptions in 
the Libor market model to that of pricing Swaptions in a multidimensional Black & Scholes (1973) lognormal model. 
Section three is then focused on finding a good pricing approximation for basket calls in this generic model. We derive 
a simple yet very precise formula where the first term is computed as the usual Black & Scholes (1973) price with an 
appropriate variance and the second term can be interpreted as approximating the expected value of the tracking error 
obtained when hedging with the approximate volatility. 

Besides its radical numerical performance compared to Monte-Carlo methods, the formula we obtain has the ad- 
vantage of expressing the price of a basket option in terms of a Black & Scholes (1973) covariance that is a linear 
form in the underlying covariance matrix. This sets the multidimensional model calibration problem as that of finding 
a positive semidefinite (covariance) matrix that satisfies a certain number of hnear constraints, in other words, the 
calibration becomes a semidefinite program. Recent advances in optimization (see Nesterov & Nemirovskii (1994) or 
Vandenberghe & Boyd (1996)) have led to algorithms which solve these problems with a complexity that is compa- 
rable to that of linear programs (see Nesterov & Todd (1998)). This means that the general multidimensional market 
covariance calibration problem can be solved very efficiently. 

The basket option representation was used in El Karoui & Lacoste (1992) where Swaptions were written as Bond 
Put options in the Linear Gauss Markov affine model. Rebonato (1998) and Rebonato (1999) detail their decom- 
position as baskets of Forwards in the Libor market model. In parallel results. Brace & Womersley (2000) used 
semidefinite programming and the order zero lognormal approximation to study the impact of the model dimension 
on Bermudan Swaptions pricing. They rely on simulation results dating back to Huynh (1994), Musiela & Rutkowski 

(1997) or lately Brace, Dun & Barton (1999) in an equity framework to justify the lognormal volatility approximation 
of the swap process and they neglect the change of measure. A big step in the same direction had also been made by 
Rebonato (1999) where the cahbration problem was reparameterized on a hypersphere. However, because it did not 
recognize the convexity of the problem, this last method could not solve the key numerical issue. In recent works. 
Singleton & Umantsev (2001) studied the effect of zero-coupon dynamics degeneracy on Swaption pricing in an affine 
term structure model while Ju (2002) use a Taylor expansion of the characteristic function to derive basket and Asian 
option approximations. 

This paper is organized around three contributions: 

• In section two, we detail the basket decomposition of Swaps and recall some important results on the market 
model of interest rates. We show that the weight's volatihty and the contribution of the forward vs. swap 
martingale measure change can be neglected when pricing Swaptions in that model. 

• In section three, we justify the classical lognormal basket option pricing approximation and compute additional 
terms in the price expansion. We also study the implications in terms of hedging and the method's precision in 
practice. 

• In section four, we expUcit the general calibration problem formulation and discuss its numerical performance 
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versus the classical methods. We specifically focus on the rank issue and its implications in derivatives pricing. 
We show how the calibration result can be stabilized in the spirit of Cont (2001) to reduce hedging transaction 
costs. 

Numerical instability has a direct cost in both unnecessary hedging portfolio rebalancing and poor risk modelling. 
By reducing the amount of numerical noise in the daily recalibration process and improving the reliability of risk- 
management computations, we hope these methods will significantly reduce hedging costs. 



2 Interest rate market dynamics 

2.1 Zero coupon bonds and the absence of arbitrage 

We begin here by quickly recalling the construction of the Libor Market Model along the lines of Brace et al. (1997). 
We note B{t,T) the discount factors (or Zero Coupon bonds) which represent the price in t of one euro paid at time T. 
We note jSj- the value at time T of one euro invested in the savings account at t (today) and continuously compounded 

with rate r^. We have pr ~ exp (^ fj^ r^ds^ . As in Heath, Jarrow & Morton (1992), to preclude arbitrage between j3r 
and an investment in the Z.C. we impose: 



B{tJ)^E^ 



T 

exp ( — / rgds 



(1) 



for some measure Q. In what follows, we will use the Musiela parametrization of the Heath et al. (1992) setup and 
the fundamental rate r{t,0) will be the continuously compounded instantaneous forward rate at time t, with duration 
9. We suppose that the zero coupon bonds follow a diffusion process driven by a t/ dimensional Q— Brownian motion 
W = {Wt,t > 0} and because of the arbitrage argument in Q, we know that the drift term of this diffusion must be 
equal to r^, hence we can write the zero coupon dynamics as: 

dB(s,T) „, 
B{s,T) 

where for all > the zero-coupon bond volatility process {<7^(f , 6); 6 > 0} is F(-adapted with values in R''. We as- 
sume that the function 9 i — > a^(f , 0) is absolutely continuous and the derivative T(f, 9) ~ d /d9{G^{t, 9)) is bounded 
on r2 X £2. All these processes are defined on the probability space {CI, {iv;f > 0},Q) where the filtration {Ft'J > 0} 
is the Q-augmentation of the natural filtration generated by the d dimensional Brownian motion W = {Wt,t > 0}. 
The absence of arbitrage condition between all zero-coupons and the savings account then amounts to impose to the 
process: 

^^''^^ -B{OJ)exp(- f a'^{sJ-s)dW,- \ f \a'^{sj ~ s)f ds] (3) 



Pt \ Jo 2 

to be a martingale under the measure Q for all T > 0. 

2.2 Libor rates, Swap rates and the Libor market model 
2.2.1 Libors and Swaps 

We note Lg [t, 9) the forwai'd 5-Libor rate, defined by: 

1 _B{t,t + 5 + 9) 
l + dLs{t,9) ~ B{t,t + d) 

and we note K{t,T) — L{t, T ~t) the forward Libor with constant maturity date (FRA). 
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A Swap rate is then defined as the fixed rate that zeroes the present value of a set of periodical exchanges of fixed 
against floating coupons on a Libor rate of given maturity at future dates T^^ and . This means: 



swap{t) 



Level {t) 



where, with cv{Ti, 7]+i) the coverage (time interval) between 7] and TJ+i computed with the appropriate basis (different 
for the floating and fixed legs) and B{t,Tf') the discount factor with maturity t/', we have defined Level{t) as the 
average of the discount factors for the fixed calendar of the Swap weighted by their associated coverage: Level{t) = 
YIi=ij cv{Tf^, T-^^ )B(t, t/^). Here t/' is the calendar for the floating leg of the swap and t/^ is the calendar for the 
fixed leg (the notation is there to highlight the fact that they don't match in general). In a representation that will be 
critically important in the pricing approximations that follow, we remark that we can write the Swaps as baskets of 
Forward Libors (see for ex. Rebonato (1998)). 

Lemma 1 We can write the Swap with floating leg Tf\ T^^ as a basket of Forwards: 

swapit) = t COi{t)K{t,T/') where aH{t) = "'^^'''lll'Jil^^'^'^'^ (4) 

with Tij. =TandO< aj,(?) < 1. 

Proof. With B{t, tP) = B{t, Tf_li){l + 5K{t,T/')), we have: 

I.UrCv{T/',T^^^)Bit,T/_l,)Kit,T/') 

swapit) = ^ — , 

Level{t) 

which is the desired representation. As the corresponding forward Libor rates are positive, we have B{t,Ti+\) < 
B{t,Ti) <B{t,Ti-i) for j e [j'r + IjA'^— 1] hence <©,(?)< 1, i.e. the weights are positive and bounded by one. ■ 

As we will see below, the weights 0)i{t) prove to have very little variance compared to their respective FRA (see 
Rebonato (1998) among others). This approximation of Swaps as baskets of Forwards with constant coefficients is the 
key factor behind the Swaption pricing methods that we detail here. 



2.2.2 The Libor market model 

As Libor rates and Swaps were gaining importance as the fundamental variables on which the market activity was 
concentrated, a set of options was created on these market rates: the Caps and Swaptions. Adapting the common 
practice taken from equity markets and the Black & Scholes (1973) framework, market operators looked for a model 
that would set the dynamics of the Libors or the Swaps as lognormal processes. Intuitively, the lognormal assumption 
on prices can be justified as the effect of a central limit theorem on returns because the prices are seen as driven by a 
sequence of independent shocks on returns. That same reasoning cannot be applied to justify the lognormality of Libor 
or Swap rates, which are rates of return themselves. The key justification behind this assumption must then probably 
be found in the legibility and familiarity of the pricing formulas that are obtained: the market quotes the options on 
Libors and Swaps in terms of their Black (1976) volatility by habit, it then naturally tries to model the dynamics of 
these rates as lognormal. 

Everything works fine when one looks at these prices and processes individually, however some major difficulties 
arise when one tries to define yield curve dynamics that jointly reproduce the lognormality of Libors and Swaps. In 
fact, it is not possible to find arbitrage free dynamics a la Heath et al. (1992) that make both Swaps and Libors lognor- 
mal under the appropriate forward measures (see Musiela & Rutkowski (1997) or Jamshidian (1997) for an extensive 
discussion of this). Here we choose to adopt the Heath et al. (1992) model structure defined in Brace et al. (1997) 
(see also Miltersen et al. (1995) or Sandmann & Sondermann (1997)) where the Libor rates are specified as lognormal 
under the appropriate forward measures but we will see in a last section that for the purpose of pricing options on 
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Swaps, one can in fact approximate the swap by a lognormal diffusion. Hence in a very reassuring conclusion on the 
model, observed empirically in Brace et al. (1999), we notice that it is in fact possible to specify Heath et al. (1992) 
dynamics that are reasonably close to the market practice, i.e. lognormal on Forwards and close (in a sense that will 
be made clear later) to lognormal on Swaps. In particular, we verify that the key property behind this approximation, 
namely the stability of the weights (Oi{t), is indeed accurately reproduced by the Libor market model. 

The model starts from the key assumption that for a given maturity 5 (for ex. 3 months) the associated forward 
Libor rate process has a log-normal volatility structure: 

dL{t,e) = {...)dt+L{t,e)Y{t,e)dW, (5) 

where the deterministic function 7 : 1 — > R'^ is bounded by some 7 e R+ and piecewise continuous. As for all 
Heath et al. (1992) based models, these dynamics are fully specified by the definition of the volatility structure and 
the forward curve today. With that in mind, we derive the appropriate zero-coupon volatility expression. Using the Ito 
formula combined with (jSj we get as in Brace et al. (1997): 

dL{t, 9) ^ {^^^ + il±^|^(T«(r, + 5){a\t^ + 5) - a\t, 0))) dt 

+ 4(1 + 5L(f , 0)) (c7^(f , + 5) - C7^(r, d))dWt 
o 

Then to get the right volatility structure we have to impose in 

a«(f , 9 + 5)- a\t, 0) = J^ML.y(t, 0) (6) 

The Libor process becomes: 

dL{t, 0) = (yj^L{t, 0) + 7(r, 0)(J^(r, + 5)L{t, 0)^ dt+L{t, 0)r{t, 0)dW, 

As in Musiela & Rutkowski (1997), we set C7^(f, 0) = for all G [0,5[ and we get, together with the recurrence 
relation ^ and for > 5 : 

With the volatility of the zero coupon defined above and the value of the forward curve today, we have fully specified 
the yield curve dynamics. 



2.3 Interest rate options: Caps and Swaptions 
2.3.1 Caps 

Let us note again j3(f), the value of the savings account. In a forward Cap on principal 1 settled in arrears at times Tj, 
j = the cash-flows are {L{Tj^i,Q) ~K)'^5 paid at time 7}. The price of the Cap at time t is then computed as: 



Capt = Y^E^ 



A(L(r,_,,0)-fc)+5 

PTj 



2.3.2 Swaptions 

To simplify the notations, we will consider that the calendars described above for the floating and the fixed legs of the 
swap are set by t/' — i5 and t/'^ = ib5, in the common case where the fixed coverage is a multiple of the floating 
coverage (for ex. quarterly floating leg, annual fixed leg). For simplicity, we will note the coverage function for the 
fixed leg of the swap as a function of the floating dates, allowing the floating dates to be used as reference in the 
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entire swap definition. From now on (7;),G[i,Af] — {Tf ),G[i,Af] and we define the coverage function for the fixed leg as 
c,5 — l{,modii=o}^^- We set ij = [5^^rJ . Using these simplified notations the Swap in @ becomes: 



swap(t) = ^(^^) B{t JN+x) ^.^j^ LeveUt) = f 5c,B(tJ,+i) 
Levelit) -tfj. 

The price of a payer Swaption with maturity T and strike k, written on this swap is then given at time t <T by: 



Swaptiorif = Ef 



L ^^^'''■^ [swapiT) - ky 



+1} 



(8) 



The expression above computes the price of the Swaption as the sum of the corresponding Swaplet prices. Because 
a Caplet is an option on a one period Swap, Caplet and Swaption prices can be computed in the same fashion. In the 
two sections that follow, we show how to rewrite this pricing expression to describe the Swaption (and the Caplet) as 
a basket option. 



2.4 Caps and Swaptions in the Libor market model 
2.4.1 Caps and the forward martingale measure 

With the Cap price computed as: 

Cap, = f B{tJj)Ej' [(L(r,_i,o) 5 

where E^i is the expectation under the forward martingale measure Qj-y defined by: 

'"^^^ = [B{QJ)^t]-' =eT{o\-Jj--)) 



dQ 

where we have noted e7'( ) the exponential martingale defined by: 

er ((T^(-, r,- - •)) = exp (^J^' CJ^(., Tj - s)dW., - ^ j^' || CJ^(., Tj -s)fds 

Let us now define the forward Libor process (or FRA) dynamics, the underlying K{t^ T) = L{t,T — t) of the Caplet 
paid at time T + which is given in the Libor market model setup in (|5} by: 

dK{tJ)^Y{t,T~t)K{t,T) [(7^(f,r-f + 5)A + rfW,] 

or again: 

dK{t, T)^y{t,T~ t)K{t, T)dW^+^ (9) 

hence K{t,T) is lognormally distributed under Pt+s- Here and in what follows, we note Vtj the cumulative variance 
from f to r and the pricing of Caplets can be done using the Black (1976) formula with Vtj equal to: 

V,j = \\y{s,T-s)fds 

Let us note that the Caplet variance used in the Black (1976) pricing formula is a linear form in the covariance. 
Recovering the same kind of result in the Swaption pricing approximation will be the key to the calibration algorithm 
design. 
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2.4.2 Swaptions and the forward swap martingale measure 

In (|8} the price of a payer Swaption is computed as the sum of the corresponding Swaplet prices, which is not the 
most appropriate format for pricing purposes. Using a change of equivalent probability measure, we now find another 
expression that is more suitable for our analysis. As in Musiela & Rutkowski (1997), we can define the forward swap 
martingale probability measure equivalent to Q^, with: 

^Q^ LUctfiiT)/liiT,+ ,) 



This equivalent probability measure corresponds to the choice of the ratio of the level payment over the savings account 
as a numeraire and the above relative bond prices are Q^— local martingale. The change of measure is identified with 
an exponential (local) Q^— martingale and we define the process ht such that: 

eT^ih.)^B{t,T)l5{T)- 



Level{t) 
which imposes: 



I=IT 



and because the volatility is bounded, we verify that e^^ (/i.) is in fact a martingale. Again as in Musiela & Rutkowski 
(1997) we can apply Girsanov's theorem to show that the process: 

is a Q'^-Brownian motion. 

Lemma 2 We can rewrite the Swaption price as: 

Swaptiont = Level{t)E^'^ [{swap{T) - (12) 

where the swap rate is a martingale under the new probability measure Q^. 

Proof. The pricing formula is a direct consequence of the change of measure above and because the Swap is 
defined by the ratio of a difference of zero-coupon prices over the level payment, it is a (local) martingale under the 
new probability measure Q'^ (below, we will see that the swap rate is in fact a Q'^— martingale). ■ 

This change of measure first detailed by Jamshidian (1997), allows to price Swaptions as classical Call options on 
a swap, under an appropriate measure. 



2.5 Swap dynamics 

We now study the dynamics of the swap rate under the Q'^ probability, looking first for an appropriate representation 
of the volatility function using the "basket of forwards" decomposition swap{t) — L/L,v <^iit)K{t, Ti) detailed in Q. 

Lemma 3 The weights (Ok{s) in the swap decomposition follow: 

N 

dcokis) = cokis) (Oiis) {a'' is, n+i -s)- a^is, Ti+y - s)) dwf 

i=ij 
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Proof. As the ratio of a zero coupon bond on the level payment and by construction of Q'^, the weights (Oi{t) must 
be Q'^— martingales (they are positive bounded). Using the forward zero-coupon dynamics, we then get: 

BltJk) ^ 5ciB(t,Ti+i) „, , r 

Level{t) .f-(^ Level{t) ^ ' '+i ' ^ 

where Wj is a Q^-Brownian motion. ■ 

We then use this result to decompose the Swap volatility as the sum of the weights volatility term and a term that 
mimics a basket volatility (the volatility of a basket with constant coefficients). We write the swap volatihty as; 

N 

dswap{s) = co,{s)K{s,Ti){Y{s,Ti-s) + ri{s,Ti))dWf (13) 

i=ij 

where the basket volatility term and the weight's residual contribution are given by: 

N f ^ \ 

£ C0iis)Kis, Ti)y{s, Ti - s) and 77 Ti) = a\s, T,+i - ^) - E ®./ W^^^^C^' Tj+i - s) 

i='T \ i='T j 

Again, the empirical stability of the weights a), (r) is the key fact at the origin of the Swaption pricing approximations 
that will follow and one of our goals below will be to show that this stability is accurately reproduced by the model. 

2.6 The forward Libors under the forward Swap measure 

We study here the dynamics of the forward Libors under the forward Swap measure. For purely technical purposes, 
we start by bounding under Q'^ the variance of the forward rates T^.), this will allow us to bound the contribution 
of the weights to the total swap variance. 

Lemma 4 With m > 1, we can bound the norm ofK{u, Tj^) under Q'^ by: 

E[K{s,ny"]<K{t,n)'"M"^(s) (u) 

where M„,[s) =exp((i — f) {inf- /2 + mf-5{N — ij)))- 
Proof. Using M \\ we can write: 

K{s, Tk) = K{t,Tk)exp /(m, 7i - u)dW^ + J' a{u, n)Y{u, n ~ u)du 

where 

N / i \ * 

a{sjk) = L <pj{s)y{sjj~s) + Y mji^j.-^) 

i=iT \./='T / '='T 

with 0,(f ) = 5K{s, 7])/(1 + 5K{s, Ti)). The corresponding forward Libor rates are positive and we have < 0,(f ) < 1 
and as in Brace et al. (1997) remark 2.3, we can bound the Forwards by a lognormal process: 

K{s, Tk) < K{t, 71) exp /(m, T^. - u)dW^, + a{u, Tk)du^ for s e [t, T] 
where we can use a convexity inequality on the norm || .||^ to obtain: 



N / i N 

L «,(^) E ^jis)7is,Tj-s)Yis,n-s) 

i=iT \j=iT J 



2 

<?>'^{l^-iTft 
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because||£jL,y (/),(f )7(5, Ti — s)y{s, — || < — hY'f, hence a{s, Tj^) = 5{N — iT)y^ which shows the desired 
result. ■ 

We now use this bound to study the impact of the weights a), (f) in the swap volatihty decomposition. 



2.7 Swaps as baskets of forwards 

For simpHcity, in what follows we will suppose that t/' = t/'^ and hence b = 1. The Swaption pricing formula that 
will be derived in section Q relies on two fundamental approximations: 

• The weights (Oi{s) for s £ [t, T] (which are Q'^-martingales) will be approximated by their value today 0), (f). 

• We will neglect the change of measure between the forward martingale measures to Q^'v+i and the forward 
Swap martingale measure Q'^. 

In this section, we study the impact of these approximations and try to quantify the pricing error they induce. The 
consequences of the first approximation are studied in lemma (|5Jl, while proposition (jS) describes the impact of the 
second. The low weight volatility is steadily observed in practice, besides Rebonato (1998), this has been studied by 
Hamy (1999) of which we report here, with the author's permission, a sample of summary statistics. The table below 
details the vol{FRA) / vol [weights) ratio in various markets, computed using the standard quadratic variation estimator 
with exponentially decaying weights (market data courtesy of BNP-Paribas London): 



Currency 

swap 
Min ratio 



USD USD GBP GBP EUR EUR 

2Y 5Y 2Y 5Y 2Y 5Y 

112 842 885 981 148 333 

Max ratio 7629 7927 6575 3473 5006 4322 

Variance .023 .020 .017 .007 .005 .004 



Sample ratio of volatility between weights and corresponding forwards. 

Here Min ratio and Max ratio are the minimum (resp. maximum) volatility ratio among the weights of a particular 
swap. We see that in this sample, the volatility of the weights is always several orders of magnitude lower than the 
volatility of the corresponding forward. Also, the weights in Q are positive, monotone and sum to one, cancelling the 
first order error terms, hence if the Forward rate curve is flat (K{s, Ti) ~ K{s, Tj) for i,j = ij, ■■■,N) we have: 



N N I 

E (0,{s)K{s,Ti)n{sJ,)^K{sJi) E (0,{s) a^isji 

i=ij i=ij \ 



s)- E (Oj{s)a\sJj-s) 

J='T 







In light of this, we will study the size of the weights' contribution to the swap volatility in terms of the slope of the 
Forward rate curve within the maturity range of the swap's floating leg. In particular, we can write the weight's part in 
the Swap's volatility as: 



E COiis)Kis,T,)iiis,T,] 

i=ij 



E (Oi{s) {K{s, Ti) - swap{s))ri{s, Ti] 



(15) 



which sets the weight's contribution as the average product of a difference of Forwards with a difference of ZC 
bond volatilities and we can expect this later term to be negligible relative to the basket volatility term in il3\ . in 
accordance with the empirical evidence. Because the payoff of the Call options under consideration are Lipschitz, we 
will approximate the Swap and forward Libor dynamics in L? under the Q"^ swap martingale measure. 
We now detail some basic properties of the weights COi{s). We note 

Lemma 5 The weights (Oi{s) defined in (0 are bounded above with: 

1 



[II •II"])'''", the L" norm. 



(Oi{s) < 



N-ij 



5swap{s) 



and satisfy \\(Oi{s)\\^ < (Oi{t) for s G [t,T]. 
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Proof. Because the weights (Oi{s) satisfy LiL,v = 1> < fi),(f) < 1 and are decreasing with / because the 
Forward rates are always positive. With: 

|a)y(i) — a),(i)| < 5swap{s) for/,y = iV,...,A^ 



we get: 



(Oiis) 



< 



1 



-5swap{s) for.se [t,T] 



and ||a),(i)||„ < ||a),(i)|| j = (Oi{t), for i e [t,T] and n > 1, because the weights are positive Q'^— martingales. | 
The next result provides a bound on the variance contribution of the weights inside the Swap rate volatility. 
Lemma 6 The L? norm of the weight's contribution in the swap volatility M3\ is bounded by: 

N 

£ C0iis)Kis,Ti)riis,Ti) 



(16) 



< 



max 1 1 (Kis, Tj ) - swa /? (i) ) 1 1 g M4 Y^S^maxj^ j^;]K{t,Tj)^{N - ir ] 



Proof. Let us note again swap{s) = YJiLij- (Oi{s)K{s, Ti), the swap rate, which we see here as the average level of 
the Forward rate curve between T and Tn- The squared norm of the weights' contribution is bounded above by: 



£ (0,{s) II {K{s, Ti) - swapis)) 77 {s, Ti) 



using a convexity inequality with YliLjj (Oi{t) = 1, < Ci),(f) < 1. To bound 77(5,74.) in this expression, we use the 
definition of G^{s, Tj^ — i) in Q and the fact that the Forwards K{s, Tj) are always positive to get: 



Y^(0i{s)[Y^5K{sJj)Y{sJj^s) 



\J=' 



with the convention Y!j=i — ~ Y!j=k if ' > ^- we recall that 7(5, Ti^ — s): R'J_ is a bounded input parameter with 

Y{s, — 5) IP < 7^, we can use (I14> and the previous lemma to get: 

I? [s, < Mtfd^naxj^^f^^^^Kit, Tj)\n - ir)' 

With these bounds we can rewrite the original inequality, using two successive Cauchy inequalities: 



^ (0,{s) {K{s,Ti)~swap{s)f\\ii{s,Ti)\\ 



< ^ ||a),(*)|l4 11(^(^,7^-) -*vvfl;,(^))|l2||77(,,7;-)||^ 

i=ij 

< max II {K{s, Tj) - swap{s))f^Mlf5^maXj^ij^^N}K{t, Tj)^{N - ij)^ 
Which gives the desired result. ■ 

With 5K{t,Tk) ~ 10^^ and {K{s, Ti) — swap{sj)~ ~ lO^-' in practice, we notice that the contribution of the weights 
to the swap volatility is several orders of magnitude below that of the basket and we will neglect it in the Swaptions 
pricing approximations that follow. Before detailing the key approximation result, we introduce some new notations. 
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Notation 7 'We define K^{s, Ti) such that: 

dK\s, Ti) = K^{s, Ti)y{s, Tf - s)dWf 
with K^{t, Ti) — K(t, Ti). We also define the following residual volatilities: 

^,{s)=K\s,Tk)y{sJk-s)-Y\s) 

with r{s) = It,^ C0i{t)K\s, T,)y{s, n - s). 

We now approximate the Swap rate with a basket of lognormal martingales parameterized by the Forward rate 
volatilities 7(5, 7]. — s) and their initial value K{t, Ti), the weights in this decomposition being equal to fi),(f ). 

Proposition 8 We can replace the Swap process by a basket of lognormal martingales weighted by constant coeffi- 
cients, with: 

^ 2' 

sup {swap{s) — Ys 

t<s<T 

< 3 niax II^X^) II4 + HK'{tJk) (N-ir) 5f)^cxp {2{T - 1) {5f {N-ir) + f/2)) 

jeliT,N] 

|2 ,^2,-2 c2.. 



+3 max \\{K{s,Tj) - swap{s))\\^Ml fd^maxj^ lirM^i' ^Tjfi^- h ) 

je[iT,N] 



where 



N 

S 



dYs^ 0}i{t)K^{sJi)y{s,T~s)dW; 

i=ij 

with Yt — swap{t). 

Proof. With the swap rate dynamics computed as in ( fT3t . we get: 

N 

d{swap{s)-Y,) - Y.^(^ki')~^k{t))K\s,Tu)y{sJk-s)dWf 

k=ij 
N 

+ £ (Ok{s) {K{s, n) - K^is, T,)) y{s, n - s)dWf 

N 

+ £ (Ok{s)Kis,Tk)riis,n)dW,^ 

We can bound the norm of the last term in this decomposition using the result in (I16t . If we look at the first term and 
note Ak,s = K{s, Tk) — K^{s, T^) with A^j = we have: 

dAk.s = A^,, (y co,{s) n -s)~ T ~ s)) ris, n - s)\ 



+K^{s, n) 1^ £ C0i{s) {a^is, n - s) - a^is, T, - s)) y{s, n - s) j ds 

+Ai,^,r{s,n-s)dw,^ 

hence: 



Akj = ^^(r, n) ^ (^M-t.^exp ^ l^k,udu^ ds^ 
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where 

N 

i=ij 

With ||/i/:.i ||2 ^ (N — It) Sf^ we can bound the norm of A/^j by: 

||A,.r II2 < K{t,n) {N-ir) 5fexp {{T - 1) {5f {N-ij) + f /l)) 

Focusing on the second term, as in ( fTsl with this time Y!i=ij O^iis) — C0i{t)—0 and ^^{s) —K^{s,Ti^)Y{s,Tf; — s) — 7"'(s), 
we can write: 

N 

^ (Wkis) - 0)k{t))K^is, Tk)7{s, n - s) 

k=ij 

The bound obtained is a function of the norm of the residual volatiHties ||^,(.?)||4 and of the spread term 7;) — swap{s))\ 

We conclude using Doob's inequality. ■ 

The term is equivalent to the variance contribution of the second factor of the covariance matrix and 

II {K{s, Ti) — sH'a/?(s))||g is a spread of rates, so we neglect both terms relative to the central volatility 7" and we 
consider the Swaption as an option on the basket F,. We notice that because we approximate one martingale by 
another, the error is in fact uniformly bounded in L^. Because of these properties and the fact that the option's payoff 
is Lipschitz, in the Swaption price approximations that follow, we will be treating the Swaption as an option on a 
basket of lognormal Forwards. 

3 Basket price approximation 

Basket options, i.e. options on a basket of goods, have become a pervasive instrument in financial engineering. Besides 
the Swaptions described in the previous section, this class of instruments includes index options and exchange options 
in the equity markets, or yield curve options and spread options in fixed income markets. In these markets, baskets 
provide raw information about the correlation between instruments which is central to the pricing of exotic derivatives. 
In this section, we detail an efficient pricing approximation technique that leads to very natural closed-form basket 
pricing formulas with excellent precision results. 

The classical "noise addition in decibels" order zero lognormal approximation was studied by Huynh (1994), 
Musiela & Rutkowski (1997) and Brace et al. (1999) when the underlying instruments follow a Black & Scholes 
(1973) like lognormal diffusion. Here, we approximate the price of a basket using stochastic expansion techniques 
similar to those used by Fournie, Lebuchoux & Touzi (1997) or Fouque, Papanicolaou & Sircar (2000) on other 
stochastic volatility problems. This provides a theoretical justification for the classical price approximation and allows 
us to compute additional terms, better accounting for the stochastic nature of the basket volatility. In fact, the first 
correction can be interpreted as a first order approximation of the hedging tracking error as defined in El Karoui, 
Jeanblanc-Picque & Shreve (1998) . 

3.1 Generic multivariate lognormal model 

We suppose that the market is composed of n risky assets S\.i— plus one riskless asset M, . We assume that 

these processes are defined on a probability space (i2,F,Q) and are adapted to the natural filtration {/v,0 < f < 7}. 
We suppose that there exists a forward martingale measure Q as defined in El Karoui, Geman & Rochet (1995) (the 
notation Q is left voluntarily non specific for our purposes here because it can either be associated with the forward 
market of maturity 7 and constructed by taking the savings account as a numeraire or it could be the level payment 
induced martingale measure as in the Swaption pricing formulas treated in the first section). In this market, the 
dynamics of the forwards are given by dF^ = F^a'^dWs and = 1 for i e [f,7], where W, is a tZ-dimensional 
Q-Brownian motion adapted to the filtration {Ft} and CTs ~ (c7^'). j ^ e jjnxJ volatility matrix and we note 

Fs e R"^" the corresponding covariance matrix defined as (Fs), ^ =< (t',<7/ >. We study the pricing of an option 



< max ll^jWllj 

^ jG[(r,A'] 
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on a basket of forwards given by Ff = Y!i=i C^iPt where 03 = (t<)()(=i n ^ '^^^ terminal payoff of this option at 
maturity T is computed as: 

for a strike price k. The key observation at the origin of the following approximations is that the basket process 
dynamics are close to lognormal. The simple formula for basket prices that we will get is specifically centered around 
a deterministic approximation of the basket volatihty: 

dFf = ( £ mi^soi] dW, with S,-,, = , (17) 

3.2 Diffusion approximation 

The classical order zero formula approximates the sum of lognormals as a lognormal variable while matching the 
two first moments. This method has its origin in the electrical engineering literature as a classic problem in signal 
processing where it represents, for example, the addition of noise in decibels (see Schwartz & Yeh (1981) among 
others). The same method was then used in finance by Huynh (1994), Musiela & Rutkowski (1997) for equity baskets 
or Brace et al. (1999) for Swaptions. Here we justify this empirical result and look for an extra term that better accounts 
for the (mildly) stochastic nature of the basket volatihty and improves the pricing approximation outside of the money. 
The approximation above simply expresses the fact that if all the forward volatility vectors were equal then the basket 
diffusion would then be exactly lognormal. It is then quite natural to look for an extra term by developing the above 
approximation around the central first-order volatihty vector JJj^ j S, ^ a/ . As in the previous section, we first define the 

residual volatihty as the difference between the original volatility and the central basket volatility Y."j=i ^j,t<^s 

and we set §j = - Y!j=\ oij,tai, for i = 1, ...,n and s e [t,T]. We also note (7® = Y!j=i ^i,tC!s (notice that is 

Ft — measurable). 

We can write the dynamics of the basket F^ in terms of a)^ and the residual volatilities ^j. Remember that for 

s G we have Qj^s > with Y!j=i = 1> hence af is a convex combination of the ai and Y!j=i ^j,s^s is a 

convex combination of the residual volatihties ^/ with Y!j=i ^j^^t = 0. As this term tends to be very small, we will 

now compute the small noise expansion of the basket Call price around such small values of Yfj=i ^j,s^s ■ We first 
write 

dFr = Fs"'' (af + eLU SijM) dWs 

(18) 

^ JS^ = ©5 (^i - erj=i ^jAs') {dWs + afds + erj=i ^jA'ds) 

and develop around small values of e > 0. As in Foumie et al. (1997), we want to evaluate the price and develop its 
series expansion in e aroimd 0. 

[F^'' -k)^\ (F,'», 6it)] with C^ = C^ + di)e + ^ + o(e2) 

We can now get the order zero term as the classical basket approximation, which corresponds to that in Huynh (1994), 
Musiela & Rutkowski (1997) or Brace & Womersley (2000). 

Proposition 9 The first term CP is given by the Black & Scholes (1973) formula. In this simple approximation, the 
basket call price is given by: 

C^=BS{T,Ff',V,,T) = F,'"N{h{Vt,T)) - KN{h{Vt,T) - ^/V;^) (19) 

where 

(ln(?) + Kr) 



h(ytj)= ^ ^ L- ^ and Vt,T= II C7, 



''^^^ds 
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where the variance can also be computed as Vtj = Tr{Q,,Ts)ds with Q,t = ojtSiJ . 

Proof. Because for s e [t,T] we have S)js > with Ly=i ^j.s = 1. as in Fournie et al. (1997) or Fouque et al. 
(2000) we can compute by solving the limit P.D.E.: 



hence the above result. Finally Tr{D.fTs) = £JLi Ly=i S}i,tO}jj < cr/, cr,' >= || erf |p allows us to rewrite the variance 
as the inner product of £2, and Fs- ■ 



We have recovered the classical order zero approximation, we can now look for an extra term by solving for C*^ 



1) 



Lemma 10 Suppose that the underlying dynamics are described by il8l . The first order term C^^\s^x,y) can be 
computed by solving: 







(20) 



j=i 2 dyj ^fi ,,=1 

= d^^fors = T 
wi thC^ =BS{s,x,Vs) given by the Black & Scholes (1973) formula as in | |79I ). 

Proof. Let us first detail expUcitly the P.D.E. followed by the price process. With the dynamics given by: 

dF,"'' = Ff'' (cjf + erj=i o^jM) dWs 

d% = (^i - el"=i Sj;,^/) (dWs + ofds + £rj=x QjAUs) 
as in Karatzas & Shreve (1991) we get for 



(Fr-^)+|(F,",a),) 



the corresponding P.D.E. : 



Lf,C' = 



= (x- yt)+ for s = T 



where is given by (with x and yt associated to F^'^ '^ and S)i,s respectively): 

2 



ds 



x^dh^ 

2 dx^ 



<.-=l 

2 



dxdyj 



xyj 



k=l 



2 dy] 



-Y^\{^i,af) + £l^y,{^^:-af,^^ 



k=l 



yr 



dyj 
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as in Fournie et al. (1997) we can differentiate this P.D.E. with respect to e to get: 




k=l 



j=i \k=i k=i 

j=l \ k=l ^ k=l 



j=i \k=i ^ 



k=l 



2\ 






j 2 dx^ 








y 



















= d'^''for5 = r 



and again as in Fournie et al. (1997) or Fouque et al. (2000) we take the limit as £ ^ oo and compute C*^'^ as the 
solution to: 

f L0cW + (L«.,3',(^/,af))x2^=0 
1 C''=0for5 = r 



which is again, with C° = BS{T,F/^ ,Vtj) given by (|T9j 

= 



or 



d') fori = r 







which is the desired result. ■ 

We can now compute a closed-form solution to the equation verified by C*^^^ using its Feynman-Kac representation. 



Proposition 11 Suppose that the underlying dynamics are described by M8\ . 
The derivative C^^^ (^t,Fj'^, icOj,t)j^i „^ can be computed as: 



/Vr.T 



(21) 



■ln^ + j;U,i,a,f)du + ^V,j 
n I ^ — : I ds 



Proof. The limiting diffusions are given by: 



F 



\(y,ffdu 



al = %,exp /■<^i^w„+/' (4/,0-oll4l K« 
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and because C^^^ solves the P.D.E. (I20> in the above lemma, with 



a^rO n(h{x,VsT)) , , ^ 1 / 1 2 

^—f = ^ ^ where n(x) = ^= exp — 



We can write the Feynman-Kac representation of the solution to ( I20t with terminal condition zero as: 



d') = [ E 



;=1 VV 



ds 



where 



M«,v) = Mll±M withy.,.= Tl 



Hence we can directly compute C^^' as: 



rT n 



'( 2 

exp(/; (cT- + ^jJ^w„j _ ^ ing + //<^w„-iy,, + iy.,/ 



ds 



which is, using the Cameron-Martin formula: 



^, « ^ (^/,(Tf)exp(2/;7^J,(T«)«f« 



Vs.T 

' In ^ + /; a^dW, + /; ( , a« ) du + \v,j 



and because for g — N{a,b 



.2\ . 



we get: 



exp 2/ (^j^O^/M 



ds 



which is the desired result. 



It is possible to compute the order two term explicitly but the computations are a bit longer We will show below 
that this result can be interpreted as a correction accounting for the misspecification of the volatility induced. 



3.3 Robustness interpretation 

The basket dynamics are essentially that of an almost lognormal process with a mildly stochastic volatility. By approx- 
imating these dynamics with a true lognormal process, we will make a small "tracking error" in hedging by computing 
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the delta using an incorrect specification of the volatility. As in El Karoui et al. (1998), we can compute this tracking 
error almost explicitly. Suppose that n<jm,.s is the value at time s of a self-financing delta hedging portfolio computed 
using the approximate volatility af . As the volatility in this delta computation is only approximately equal to the 
volatility driving the underlying assets, there will be a small hedging tracking error Cs computed as — 4 — n^o) 5 
for s e [t,T], where P^f is the price of the option at time s, computed using the approximate volatility (xf . Of course, 
we know that Pa^,T — (F/" — K)^ and we can understand E[es] as a price correction accounting for the volatility 
misspecification. From El Karoui et al. (1998) we know that we can compute this (exact) tracking error explicitly as: 













( 


n 


-\K\n 












!=1 





From the computation of C''' in the previous part we know: 

rT 



ds 



C 



(1) - 



.7=1 



, n{h{V,,,T,Fn) 



ds 



With o\ = af + £,^, and because ^"^1 S,,, = 1, we rewrite (I22t as: 



2Jt 




,,2 d^C\F-,V,j) 



d^C\F^,V,j) 



ds 



ds 



The first order expansion of ej for small values of gives: 



'T 



ds 



(22) 



writing the value of the Gamma expHcitly, we get: 



J n 



Ah{Vs,T,Ff)) 



ds 



and finally C^^^ ~ E ej ^ . This means that the first order correction in the basket price approximation can also 
be interpreted as the expected value of the first order tracking error approximation for small values of the residual 
volatility (^^'. This validates the price approximation in terms of both pricing and hedging performance. To make the 
link with section two explicit, we now write the order zero approximation in the particular case of Swaption pricing. 



3.4 Swaption price approximation 

If we go back to the particular Swaption pricing problem developed in section two, the result above allows us to 
approximate the price of a Swaption. 

Proposition 12 Using the above approximations, the price of a payer Swaption with maturity T and strike K, written 
on a Forward Swap starting at T with maturity Tn is given at time t <T by the Black formula plus a correction term: 



Swaption, = Level{t) {swap(t)N{h) — KN{h — ^Vtj)) +Level{t)C 



1) 



(23) 



with 



h = 



In 



swap{T) 



with Vt,T = 



\ris)\\'ds 
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where swap(t, T, 7}v) is the market value of the Forward swap today with 



and 



d),{t) = ©KO-^^ and fis) - f a},{t)r{s,Ti-s) 

swapit) 



CW = f t %(0 ^^"^^^^ exp ( 2 / {^is,Tj^s),f>is))du 

Jt \/Vt.T 




'\ni^+j:{^{s,T,-s),f>{s))du + ^^V, 



where t, [s, Ti - s) = 7(5, Tj-s)- y^{s). 

Proof. This is a direct consequence of proposition and (II It . with F!. — K^{s, Ti) and = 7(5, 7) — s). ■ 

In the last section, we will study the practical precision of this approximation by comparing the price obtained 
using the formulas above with the price obtained by Monte-Carlo simulations in both the Libor Market model and in 
the generic multidimensional Black & Scholes (1973) model. 

4 Libor market model calibration 

In this section, we detail the calibration problem and its resolution by semidefinite programming techniques. For a 
general overview of semidefinite programming algorithms see Nesterov & Nemirovskii (1994) or Vandenberghe & 
Boyd (1996). Because it provides sufficient precision in most market conditions, we will use the order zero approxi- 
mation here (if the rates become less correlated and the relative variance of the second factor increases, we can always 
replace i2, below by a new matrix, factoring in the first order price correction). Let us write the market variance in the 
approximation obtained in the last section as a function of the scalar product of the forward rates covariance matrix 
and a matrix computed from market data on the Swap weights: 



V,,T = 



N 

Y,a)i{t)Y{s,T~s) 



fT ( N N \ 

= / ( E E (7(^, Ti - s),y{s, Tj ~ s))\ ds 



Tr{D.,Xs)ds (24) 

where f2,,X, e R^^^,f e [0, T] are positive semidefinite symmetric matrixes defined by: 

n, = cb{t)(bitf - iMt)&jit)),,je[i,N] ^ and X, = {{y{s,Ti - s),Yi^,Tj - s))) ^.^^^j^^ >z 

i.e. Xs is the covariance matrix of the forward rates (Gram matrix of the 7(5, 7] — s) vectors). This shows that the cu- 
mulative market variance of a particular Swaption can be written as a linear functional of the Forward rates covariance 
matrix. With (y^j^^f^g., T/^ for k— 1 , . . . , M, the market cumulative variance for the Swaption of maturity 7^. as inputs, the 
calibration problem can then be written as an infinite-dimensional linear matrix inequality (L.M.I.) : 

Find X, 



s.t. Tr[ai (^// X,ds) j - (t2^,,,,,,71 for^= (1,...,M) (25) 
XshO 

in the variable X^ : R+ S„, where the matrix is quoted by the market today. 

Because the market variance constraints are Unear with respect to the underlying variable X^ and the set of positive 
semidefinite matrixes is a convex cone, we find that the general caUbration problem is convex and given a convex 
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objective function, it has a unique global solution. For simplicity now and to keep the focus on the problem geometry, 
we discretize Xs with a 8 frequency and make the common (but not necessary here) simplifying assumption that 
although the forward rates volatilities are not stationary, their instantaneous correlation is, hence the volatility function 
take a quasi- stationary form y{s,x) — a{s)r\{x) with a and 77 such that a{s) — (t(^ [^■'J): ^(") = ^(i L^"J) 
o{s) = Tj (s) = when s < 0. The expression of the market cumulative variance then becomes Vr.r — Tj=r ^Tr {Q.tXi). 
We can account for Bid-Ask spreads in the market data by relaxing the constraints as: 

Find Xi 

s.t. al, ,n < If=, 8Tr {a,^,Xi) < al., ,n for k = 1, ...,M (26) 
X;^Ofor/ = 0,...,r 

where we have set — ((7^(5) (t7(7^' — i), Tj(ry — 5))) . ^.^jj n]—^ (keeping in mind that the vectors r]{Ti~s) creating 
this matrix "shift" from period to period). Numerical packages such as SEDUMI by Sturm (1999) (for symmetric cone 
programming) solve these problems with excellent complexity bounds similar to those obtained for linear programs 
(see Nesterov & Todd (1998)). 

4.1 Applications 

In general, the calibration problem gives an entire set of solutions. Different choices of convex objectives are detailed 
below. 

4.1.1 Bounds on other Swaptions 

One of the most simple choices of objective matrix C is to set it to another Swaptions associated matrix £2?; . The 
calibration problem finds the parameters for the Libor market model that gives either a minimum or a maximum 
arbitrage-free price (within the BGM framework) to the considered Swaption while matching a certain set of market 
prices on other Caps and Swaptions (see d' Aspremont (2002^?) and d'Aspremont (2002a)). 

4.1.2 Distance to a target covariance matrix 

Let A be a target covariance matrix (for example, a previous calibration result or an historical estimate), we can 
minimize \\A —X\\ under the constraints in j26> . If ||.|| is the spectral or Euclidean norm, this is a symmetric cone 
program and can be solved as in Nesterov & Todd (1998) or Sturm (1999). 

4.1.3 Maximum entropy 

In the spirit of Avellaneda, Friedman, Holmes & Samperi (1987), let P be a covariance matrix representing prior 
information on the distribution of Forwards, as in Vandenberghe, Boyd & Wu (1998) we can minimize — Indet (X) + 
Tr{p-^X) to find the maximum relative entropy solution to the cahbration problem. 

4.1.4 Smoothness constraints 

It is sometimes desirable to impose smoothness objectives on the calibration problem to reflect the fact that market 
operators will tend to price similarly the variance of two products with close characteristics. A common way of 
smoothing the solution is to minimize the surface of the covariance matrix that we approximate here by: 

S= \\^'J^\\^ where AijX ^ 

ijell.n] 

Again, this is a symmetric cone program. 



X, 
X, 



-Xi-ij 
-Xij-i 
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4.1.5 Calibration stabilization: a Tikhonov regularization 



Along the lines of Cont (2001), we can explore the impact of the smoothness constraints introduced above. We can 
think of the calibration as an ill-posed inverse problem and write the smooth calibration program as a Tikhonov (1963) 

regularization of the original problem. If we set, CJ, 



2 

Mid.k ■ 



Bid.k ' 



■'Aks.k 



/2, minimizing 



'Mid,k'^k 



ol\\X\\ will then directly improve the stability of the calibration problem. 



4.2 Rank Minimization 

Because the calibrated model will be used to compute prices of other derivatives using mostly Monte-Carlo techniques 
or trees, it is highly desirable to get a low rank solution. In general, the matrix solution to the calibration problem will 
lie on the border of the semidefinite cone and hence will be singular but there is no guarantee that the rank will remain 
below a certain level. In general (cf.Vandenberghe & Boyd (1996)), this problem is NP-Hard. However, some very 
efficient heuristical methods (see Fazel, Hindi & Boyd (2000) on trace minimization) can produce results with very 
rapidly decreasing eigenvalues. In practice and in accordance with prior empirical studies (see Brace et al. (1997)), 
all solutions (even those with a high rank) tend to have only one or two dominant eigenvalues with the rest of the 
spectrum several orders of magnitude smaller. 



5 Numerical examples 
5.1 Approximation precision 

To assess the practical performance of the lognormal swap rate approximation in the pricing of Swaptions, we will 
compare the prices obtained for a large set of key liquid Swaptions using Monte-Carlo simulation and the lognormal 
forward swap approximation. We have used the classic Euler discretization scheme as detailed for example in Sidenius 
(1998). In figure ([0, we present a plot of the difference between two distinct sets of Swaption prices in the Libor 
Market Model. One is obtained by Monte-Carlo simulation using enough steps to make the 95% confidence margin of 
error always less than Ibp. The second set of prices is computed using the order zero approximation formula above. 
We can notice that the absolute error is increasing in the underlying maturity of the Swaption and that its sign is not 
constant. This plot is based on the prices obtained by calibrating the model to EURO Swaption prices on November 

6 2000 (data courtesy of Paribas Capital Markets, London). We have used all Cap volatilities and the following 
Swaptions: 2Y into 5Y, 5Y into 5Y, 5Y into 2Y, lOY into 5Y, 7Y into 5Y, lOY into 2Y, lOY into 7Y, 2Y into 2Y, lY 
into 9Y (the motivation behind this choice of Swaptions is liquidity, all Swaptions in the lOY diagonal or in 2Y, 5Y, 
7Y, lOY are supposed to be more liquid). The absolute error is always less than 4 bp which is significantly lower than 
the Bid-Ask spreads. 

In the second figure we plot the error in the basket pricing formula for a basket of assets, having supposed 
that the forwards are all martingale under the same probability measure (hence we test the precision of the approxi- 
mations without the error coming from the forward measures, this is also a test of the formula's precision in an equity 
framework). The reference is given by a Monte-Carlo estimate with 40000 steps. The numerical values used here are 
= {0.7,0.5,0.4,0.4,0.4}, O), = {0.2,0.2,0.2,0.2,0.2}, T ^5 years, and the covariance matrix is given by: 
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/ 0.64 0.59 0.32 0.12 0.06 \ 

0.59 1 0.67 0.28 0.13 

0.32 0.67 0.64 0.29 0.14 

0.12 0.28 0.29 0.36 0.11 

V 0.06 0.13 0.14 0.11 0.16 / 



The covariance used here comes from an historical estimate and has the typical level, spread, convexity eigenvector 
structure. These values are meant to replicate the pricing of a 5Y into 5Y Swaption without the change in measure. 
We can see that the pricing error is less than 2bp with the order zero approx. and the additional order one term does 
not provide a significant benefit. In fact, the order zero term reaches an excellent precision near the money, a feature 
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Absolute pricing error (in basis points) 



5 T- 




Swaption (Maturity, Underlying) 



Figure 1 : Absolute error in the order zero price approximation versus the Libor market model prices estimated using 
Monte-Carlo simulation, for various ATM Swaptions. 

that is constantly observed when the covariance matrix has the structure given above, where the first level eigenvector 
accounts for around 90% of the volatility and the model is close to univariate (as noted in Brace et al. (1997)). 
However, we observe in figure (|3} that the order one approximation does provide a significant precision improvement 
when the rates are less correlated. Finally, in a pure equity case, i.e. when the initial value of the underlying assets 
is not significantly smaller than one (an equity basket option for example), the order one correction very significantly 
reduces the relative error, as can be observed in figure ©. 

5.2 Calibration 

Using the same data set as above, we calibrate a covariance matrix under smoothness constraints. The resulting matrix 
is plotted in figure (|5|i. In figure ^ we plot the eigenvectors of this matrix. The first vector has a level shape while 
the second one is close to a spread of rates. We can notice that this purely market implied covariance factor structure 
closely matches the results obtained using estimates from historical data. 

6 Conclusion 

The methods described in this work are organized around one central objective: the design of a true "black-box" cali- 
bration and risk-management tool for classic multifactor interest rate models. In particular, the performance guarantee 
given by the numerical methods used here makes it possible to design a calibration procedure that does not require nu- 
merical baby-sitting. Furthermore, the possibility of stabilizing the calibration result should induce significant savings 
in hedging transaction costs by suppressing the possibility of purely numerical calibration hedging and hence P&L 
hikes. 
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Figure 2: Order zero (dashed) and order 
one (plain) absolute approximation error ver- 
sus the multidimensional Black-Scholes bas- 
ket prices obtained by simulation for various 
strikes. 



Figure 3: Order zero (dashed) and order 
one (plain) absolute approximation error ver- 
sus the multidimensional Black-Scholes bas- 
ket prices obtained by simulation for various 
strikes. (Diagonal covariance matrix) 
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zero (dashed) and order one (plain) relative 
approximation error versus the multidimen- 
sional Black-Scholes basket prices obtained Figure 5: Smooth calibrated covariance ma- 
by simulation for various strikes. trix. 
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Figure 6: First eigenvector "level" Figure 7: Second eigenvector "spread" 
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In practice however, two important obstacles remain in the design of a "Swiss army knife" interest rate model: 
smile modelhng and rank reduction. It is at this point not possible to globally cahbrate the model to both the smile and 
the covariance structure, instead, one has to apply a two-step procedure to first cahbrate the correct smile structure and 
then recover the covariance information using the methods detailed here. This makes it impossible to jointly optimize 
the cahbration result on the smile and the covariance structure (for smoothness, stability, etc.). The second problem 
is rank reduction: numerical methods for American-style securities pricing are only efficient for models with a small 
number of factors. Empirical evidence suggests that market covariance matrixes should have a rapidly decreasing 
eigenvalues and the matrix cahbrated from market data on Caps and Swaptions display that behavior, hence the rank 
reduction is essentially a numerical backward compatibility problem and recent advances in quantization methods (see 
Bally & Pages (2000)) or American Monte-Carlo (see Longstaff & Schwartz (1998) for example) make it reasonable 
to believe that this hmitation will eventually be hfted. 
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